Theory of finite-entanglement scafing at one-dimensional quantum critical points 
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Studies of entanglement in many-particle systems suggest that most quantum critical ground 
states have infinitely more entanglement than non-critical states. Standard algorithms for one- 
dimensional many-particle systems construct model states with limited entanglement, which are a 
worse approximation to quantum critical states than to others. We give a quantitative theory of 
previously observed scaling behavior resulting from finite entanglement at quantum criticality: the 
scaling theory of finite entanglement is only superficially similar to finite-size scaling, and has a 
different physical origin. We find that finite-entanglement scaling is governed not by the scaling 
dimension of an operator but by the "central charge" of the critical point, which counts its uni- 
versal degrees of freedom. An important ingredient is the recently obtained universal distribution 
of density- matrix eigenvalues at a critical point The parameter- free theory is checked against 
numerical scaling at several quantum critical points. 

PACS numbers: 64.70.Tg, 03.67.Mn,75.10.Pq 



A system in its ground state can undergo a second- 
order (continuous) quantum phase transition as a control 
parameter is varied through a critical value. As at a ther- 
mal phase transition, the critical point has correlations 
over long length scales. These correlations are described 
by properties of the critical point that are "universal" , 
i.e., independent of microscopic details. The entangle- 
ment entropy is a measure of the quantum-mechanical 
nature of correlations and in many cases is also univer- 
sal 0, i, i, i, i, S S • The entanglement entropy of a 
pure state of a bipartite system AB is defined as 



S 



-TipA^ogPA = -TrpB logpB, 



(1) 



where pA (pb) is the reduced density matrix of subsystem 
A{B). A nonzero entanglement entropy reflects that in 
quantum mechanics, a "complete description" of a sys- 
tem (i.e., a pure state) does not imply a complete de- 
scription of subsystems. 

We would like to understand a consequence of diverg- 
ing entanglement at quantum criticality: any approach 
constructing states with a limited amount of entangle- 
ment will show universal, systematic errors in describ- 
ing quantum critical states. It was shown previously 
by Tagliacozzo et al. [lol | in a numerical study of two 
one-dimensional critical models that finite entanglement 
leads to scaling behavior like that induced by other per- 
turbations of a critical point: it introduces a finite corre- 
lation length ^ ~ x**! where x (defined below) is related 
to how much entanglement is retained and k is the finite- 
entanglement scaling exponent. That work gave convinc- 
ing evidence for this behavior in the two models studied 
but did not attempt to explain its origin or develop a 
theory predicting k. As the retained entanglement x in- 
creases, ^ — > oo and criticality is restored. 

The main result of this Letter is a theory for this be- 
havior for conformally invariant critical points in one di- 
mension. The theory predicts that k, unlike other scal- 



ing exponents, is determined by the central charge of the 
critical point. It leads to a specific formula for this de- 
pendence and also explains the observed scaling of entan- 
glement entropy. The iTEBD algorithm ll|, which has 
finite-entanglement errors but not finite-size ones, is then 
used to study a number of critical points. The numeri- 
cal results confirm that central charge determines finite- 
entanglement scaling and show parameter-free agreement 
with the theoretical predictions. 

The best understood quantum critical points are in 
one spatial dimension, where most translation-invariant 
systems have critical points with "conformal invariance" , 
an infinite symmetry group related to conformal maps of 
the plane. The entanglement entropy between two halves 
of a large one- dimensional system close to the critical 
point (large correlation length ^) is [1]: 



S 



glog(?/a) 



(2) 



Here c is the "central charge" of the critical point, a 
number that counts how many degrees of freedom of the 
system are critical, and a is a short-distance length scale, 
such as the lattice spacing in a spin chain. 

Our goal is to understand how the increased entan- 
glement of a quantum critical state (Eqn. ([2])) affects 
classical representation of the state, e.g., in a computer 
algorithm. A powerful approach to obtain low-energy 
states of one-dimensional systems is the density-matrix 
renormalization group (DMRG) algorithm [l^ . This al- 
gorithm and its descendants [ll| construct trial wave- 
fimctions that are "matrix product states" (MPS) [ist . 
Consider a system with N sites, where each site has d or- 
thogonal states. Any pure state of the system is a super- 
position of the product basis states \s1S2 ■ ■ ■ sn) = \{s}) 
where 1 < < d. A MPS for such a system has the form 
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For each site i there are d matrices A*, of a finite di- 
mension X ^ X- A product wavefunction describing an 
unentangled chain of particles/spins is obtained by mul- 
tiplying together scalar amplitudes for the particle/spin 
state at each site. The MPS generates entanglement by 
using matrices instead of amplitudes for each site, and 
an increasing amount of information can be stored as the 
matrix dimension x increases. 

A convenient representation of this state \'^) for entan- 
glement purposes is the Schmidt decomposition. Split- 
ting the system into two parts at one bond, the Schmidt 
decomposition is a basis choice that expresses the original 
wavefunction as a sum of product states of wavefunctions 
for the two halves of the system: 



(4) 



The |$„a) and |<f>„B) form orthonormal bases for the 
Hilbcrt spaces of the subsystems A and B to the left 
and right of the bond. The Schmidt decomposition con- 
tains more than one term for entangled states such as 
an Einstein-Podolsky- Rosen state : measurements on 
subsystem B put the system in a state with prob- 
ability and force the left half of the system into the 
corresponding basis state. This converts the system into 
a mixed state with entropy 



S = -J2xl\ogXl 



(5) 



We use the recently developed infinite Time-Evolving 
Block Decimation (iTEBD) algorithm llj to study sev- 
eral local one-dimcnsional Hamiltonians with translation- 
ally invariant ground states. Exploiting translational in- 
variance, there are only finitely many different matrices 
Al , which are found by minimising the energy of the 
wave function in Eqn. Because this method always 
constructs wavefunctions for the infinite system, its er- 
rors result from finite entanglement rather than finite 
size; it was similarly used in [lO| for the original numer- 
ical study of how finite entanglement affects quantum 
criticality. Here our focus is a theory for how the nu- 
merical observations differ fundamentally from conven- 
tional finite-size scaling, together with additional iTEBD 
results. We study cases in which the iTEBD algorithm 
converges effectively to the matrices of a given dimension 
that minimize the expectation value of the Hamiltonian 
llj . Away from a critical point, the entanglement en- 
tropy of the exact wavefunction is finite, and the approx- 
imation by matrix product states converges rapidly [isj . 

At the critical point, finite-dimensional matrices can- 
not approximate the ground state as well; the maximum 
entropy for x Schmidt eigenvalues is logx- The numer- 
ical eigenvalues for the system at the critical point but 
with finite x actually obey a distribution that docs not 



maximize the entropy. Instead, there is a physical inter- 
pretation of the observed distribution, which underlies 
the theoretical analysis. We first note that the observed 
eigenvalue distribution at finite x and at the critical point 
(correlation length ^ = oo) is similar to the infinite-x dis- 
tribution away from the critical point (at some finite ^ 
determined by %). This is the microscopic origin of the 
observation [10| that physical quantities computed at in- 
finite ^ and finite x reproduce those at finite ^ and infinite 
X, with an empirical scaling law ^ oc x'^ in which the ex- 
ponent K depends on the specific critical point. More pre- 
cisely, the eigenvalue distribution of the reduced density 
matrix of a one-dimcnsional system with a large correla- 
tion length ^ was recently studied [l[: the mean number 
of eigenvalues larger than a given value A is 

(6) 



n(A, b) = /o(2V&(-&-21ogA)). 



Here Iq is the zeroth modified Bessel function and b is 
a parameter that varies with g, diverging at the critical 
point; in fact, summing Eqn. ([5]) shows that b is deter- 
mined by the entanglement entropy: 6 = S/2. We find 
that this distribution is a good description of the ob- 
served eigenvalue distributions, even at criticality, with 
finite but large x- Wc now discuss a theory for the finite- 
entanglement scaling exponent k, which also appears in 
the scaling of entanglement entropy with x, 



6 6 



(7) 



We wish to compute the effective correlation length ^ 
that describes the state obtained by minimization of en- 
ergy at the critical point and finite x- Since the true 
ground state has infinite ^, a state with a finite ^ car- 
ries an energy cost proportional to 1/^^ (We as- 
sume that the critical point separates two gapped phases, 
rather than terminating a critical phase.) However, this 
energy cost is balanced because the retention of only the 
first X eigenvalues has a more severe effect on a ^ = oo 
critical state than a noncritical state. Hence the energy 
(E) for finite X: as a function of correlation length ^, 
should include a second term that computes the energy 
cost from discarding all but the largest x eigenvalues: 



E{0 = Eo 



A 

e 



Pr{b,x) 



(8) 



Here Eao is the ground-state energy of the critical sys- 
tem per unit length and Pr{b, x) = Si^x ^^"^ 
residual probability, making the b dependence explicit. 
Note that (3 and A are non-universal constants that will 
drop out of the asymptotic behavior for large x [isj - 

In the asymptotic large-x limit, the discarded eigen- 
values can be assumed to form a continuum and be dis- 
tributed according to Eqn. The residual probability 
is now an integral with the values of A(6, n) obtained by 
inverting Eqn. For large-x, we find 



Pr{b,x) 



logx - 2& 



X -(i°gx)^ 
e IS 



(9) 
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Replacing ^ in favor of b using ^ = e^^^'^ = gi2f)/c ^^^^ 
using Eqn. (|31|) in Eqn. ([8]) gives the energy per unit 
length as a function of b and x- To find the ground state 
at fixed Xj we minimise the energy with respect to b and 
find that 



6 



12+1 



^5 = 



1 



/12 +1 



■logx, 



(10) 



with corrections of order 1/logx- The relationship be- 
tween K and c is a central result of this Letter. 

We have performed numerical tests of the scaling pre- 
diction in Eqn. (|35p on several critical points. The spe- 
cific one-dimensional systems we have studied are the 
quantum Ising model in a transverse field, the XX Z spin 
chain for spin-1/2, and the spin-1 generalized Heisenberg 
model. The quantum Ising model whose Hamiltonian is 



i+l 



is critical at g = 1 with central charge c = 1/2. 
XX Z spin chain is described by the Hamiltonian 
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FIG. 1: Scaling of the entropy with matrix dimension x'- The 
entropy calculated numerically using the iTEBD algorithm for 
several critical points in one dimension: the XX and Heisen- 
berg (XXX) models, the transverse Ising model, and the s — 1 
biquadratic Heisenberg model (cf. Methods). It can be seen 
that in each case S oc log x- Furthermore, the lines for models 
with the same central charge c have the same slope, although 
other universal properties differ. 



and is critical in the entire range 7 £ [0, 1] with central 
charge c = 1. Critical exponents change continuously 
with 7 [l3| ■ The last system we have studied is the S = 1 
Heisenberg chain with biquadratic term, 

H = Y, [cos0(S,.S,+i) + sin0(S,.S,+i)2] . 



This system is known to have two exactly solvable critical 
points, the SU(2)2 point at 61 = -7r/4 with c = 3/2 ^ 
and the SU(3)i point at 6* = 7r/4 with c = 2 The 



entire region S [7r/4, 7r/2) is critical with c = 2 [2lj . 
which is consistent with our results. A detailed study of 
critical S\J{N) can be found in Ref. 

The prediction that the scaling of the entropy S 
depends only on c can be checked directly (Fig. [T]). 
The relationship between the entropy and the largest 
density-matrix eigenvalue, which is used by Calabrese 
and Lefevre [if, is approximately satisfied (Fig. (2]). The 
model above implies that the energy in the critical Hamil- 
tonian for a state with finite value of ^ should go as 
E = Eq + B/S^'^, where Eq is the actual ground-state 
energy and B a non- universal constant. This connection 
between the actual correlation length ^ and the energy 
was tested (Fig. [3]) for ground states on and off criticality. 

Numerical calculations can estimate the entropy scal- 
ing using the definition of the entropy in Eqn. ([2]) as well 



23[, and compare 



as the relations Eqn. ^ and S = 2b 
the results to the parameter-free theoretical prediction 
Eqn. (|35p. The results are shown in Fig. 01 The agree- 
ment of the numerical values among themselves and with 
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FIG. 2: Scaling of the parameter b: b = — log Xmax calculated 
numerically for different models, where Xmax is the largest 
eigenvalue of the Schmidt decomposition. For all models, b oc 
logx for sufficiently large matrix dimension x- Comparison 
to Fig. [1] shows that the relationship S = 26 is well satisfied. 



Eqn. ((35)) is at worst about 20%; even for small x the er- 
ror in the asymptotic theoretical prediction is comparable 
to the discrepancies between definitions of the entropy. A 
stringent test of the nonlinear c dependence of Eqn. ([35]) 
is that systems of decoupled copies arc also found numer- 
ically to obey this scaling prediction (not shown). 

These results suggest that the effect of finite entan- 
glement in any MPS description near a quantum criti- 
cal point in one dimension results from universal proper- 
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FIG. 3: Energy vs. correlation length: The energy calculated 
numerically scales as 1/^^ with correlation length ^. The 
same scaling also appears (inset) when the state in question 
is obtained as the ground state of a slightly ofT-critical Hamil- 
tonian. The correlation length in a matrix product state can 
be obtained from the transfer operator [20| . 



tics of the quantum critical point, specifically its central 
charge rather than a scaling dimension as in finite-size 
scaling. There are several potential extensions of this ap- 
proach. In addition to studying other types of quantum 
critical points, there are recently developed MPS-type 
algorithms for time-dependent problems [25| and higher 
dimensions [2^ that may also show universal scaling of 
errors near quantum criticality. 
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FIG. 4: K for different values of the central charge: k cal- 
culated numerically using different definitions of the entropy 
(equations ^ nearly agree ((S}, while S = 2b gives different 
values as shown) , and according to the parameter-free asymp- 
totic theoretical result in Eq. (|35p . Numerical evaluations of 
the theoretical model for finite x sue also shown. One cause 
of discrepancies between the values corresponding to different 
definitions of the entropy is that only finite x can be stud- 
ied; the different definitions of entropy should converge in the 
asymptotic large-x limit. The numerical values of k at c = 1/2 
and c — 1 are consistent with those reported in Ref. [l^ 
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Supplementary information: Here we outline the 
calculation of the asymptotic scaling form of the entan- 
glement entropy. The procedure involves writing down 
the energy density (energy per unit length) of the criti- 
cal system for states which have only a finite number (%) 
of terms in their Schmidt representation. We argue that 
as long as we are in the scaling limit, all states and hence 
their corresponding energy density for a fixed x depends 
on only one other parameter for a given system, which 
can be chosen so that the energy density is minimised. 
The state thus obtained is the one for which we calculate 
the entanglement entropy to investigate its scaling with 
X- 

The entanglement entropy (S) in the vicinity of the 
critical point is given by 



fi 



6 



(11) 



where c is the central charge and ^, the correlation length. 
The other definition of the entropy is [l[ 



S = 2b, 



(12) 



where b = — 21ogA„ia2; and A^„q^ is the largest eigen- 
value of the reduced density matrix, b is also the param- 
eter that appears in the Calabrese-Lefevre distribution 
for the eigenvalues of the reduced density matrix in 
the vicinity of the critical point [l[ , 



n{X) = lo (2V-fo2-261ogA 



(13) 



Io{x) is the zeroth modified Besscl function. From 
Eqns. [Tl] and [121 we obtain a relation between b and 
^, which is 



(14) 



The leading correction to the free energy density (F) 
of a critical point with conformal invariance is (cf. I. 
Affleck, Phys. Rev. Lett. 56, 746 (1986)) 



6v ' 



(15) 



where is the energy density of the critical state and 
V, the velocity of low energy excitations. Using the ther- 
modynamic identity 



E = F-T 



dF 
df 



(16) 



we can calculate the energy density (E) at temperature 
T: 



E^E^ 



ncT' 
6v 



(17) 



For a one dimensional critical point with conformal in- 
variance, a non-zero temperature T is equivalent to a 



finite correlation length ^ with ^ cx 1/T. Thus, the en- 
ergy density {E) of a state with a finite ^ at the critical 
point of a system is given by 

E^E^ + ^, (18) 

where A is a non- universal constant. The critical state 
has infinite entanglement entropy and thus requires an 
infinite number of terms in the Schmidt decomposition 
to represent. To study the scaling of the entanglement 
entropy, we allow only a finite number (x) of terms in the 
representation of the ground state of the critical system 
and calculate the value of S' as x ^ cxd. The effect of a 
finite X is to introduce a finite correlation length ^. This 
implies that the spectrum has a gap A ^ 1/^. We now 
argue that another term needs to be added to Eqn. [181 
for states with a finite x- To see this we first note that the 
state IV'o) which has the energy density given by Eqn.fTSl 
for a given ^ has the Schmidt decomposition 



(19) 



where and \ipn) are states of the two (semi- infinite 
halves) of the system to the left and right of a given point 
respectively. A„ has a distribution like in Eqn. [T31 with b 
given in terms of ^ by Eqn. 1141 and with the normalization 



E 

n=l 



A 



(20) 



n in Eqns. [19] and [20] goes all the way to infinity as re- 
quired by Eqn. 1131 However, the state with finite x is 
a truncation of the sum in Eqn. [19] with only x terms. 
This state \ip) which is an approximation to {ipo) is thus 



\/Yln=l 



(21) 



and has an energy density different from the one given 
by Eqn. [18] A knowledge of the energy density Eq of the 
ground state |V'o) and the excitation spectrum, permits 
us to calculate the energy density E oi \ip) as 



E 



i^oKV'olV')!' 

Eq + {Eex - 



+ Eex ( 

- Eo) (1 



i(^oiV')n. 



(22) 



Eex is a measure of the energy density of the excited 
states and thus {E^x — Eq) ^ A ^ 1/^. Using this fact 
and Eqns. [T9l and [21] in Eqn.[22l we obtain 

E = Eo + (^(^l-p^Xl}j ^Eo + ^Prix). (23) 

In the above equation /3 is a non-universal constant and 
the residual probability Prix) is given by 



n=x+l 



A„ 



(24) 
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The p term indeed describes an energy cost per length 
because we have computed the energy cost per trunca- 
tion, and the number of truncations is hnear in length. 
We can combine Egns. [T51 and E51 to obtain the energy 
per length of the ground state of the critical system for 
finite x, 



(25) 



There is an optimum value of ^ which mimimizcs E for 
a given x- If £. is very large, the second term on the right 
hand side of Eqn. [55] is very small. However, this also 
means that b is very large from Eqn. [Til and the A„'s fall 
off very slowly with n from Eqn. 1131 Thus, the residual 
probability after discarding a finite number of terms is 
large and so is the third term on the right hand side of 
Eqn. [55] in comparison to the second term. On the other 
hand, if ^ is very small, the effect is the opposite with the 
second term becoming large and the third small. Thus, 
there is an optimum ^ (and hence b) for a given x that 
minimises the energy density. We will determine this 
value of b in the limit x — > cx). In this limit we can make 
n a continuous variable and the residual probability 



Prib,x) 



X{b, n)^dn, 



(26) 



We will see from the final solution that b cx logx and 
thus '°^2v^^^ indeed large when x is large. 

We can now determine which state \tjj) minimises the 
energy E in Eqn. [25] lip) will depend on the x values of 
A (6, n), which in turn depend on the single parameter b as 
does ^ through Eqn. 1141 Thus, it is sufficient to minimise 
the energy with respect to b. Taking the derivative of the 
energy in Eqn. [25] and setting it equal to zero gives 



I 2 log X 1 _ 
"^(logx-2h)^ / 



(logx) 



462 

e " 



2b 



log X— 2b 

(32) 



In the limit of large x, the above equation has a solution 
of the type 



b = Mlogx- 



(33) 



Substituting /i from the above equation in Eqn.[32]lets us 
calculate its value. The first term in the curly brackets 
is independent of x- while the second term is 0{ ^J^^ ) 
and can be neglected. Thus, the two leading terms are 
both of the form x to an exponent and the two exponents 
have to be equal for Eqn. [32] to hold. This gives us the 
following quadratic equation for /i: 



where we have made explicit the fact that the eigenvalues 
and hence also the residual probability depend on b. We 
can invert Eqn. [T5lto obtain X{b,n). If x is large, so is n 
and we can use the asymptotic form of Eqn. 1131 which is 



(1 - 2^J) 



2 48^2 



which has the solution 



(34) 



n{X) 



„2^-b2-2blogA 



47rv/-&2 - 26 log A 



(27) 



Thus 



logn(A) = 2V-&^ - 261ogA+-log (^^TTyJ-b^ - 26 log A 

(28) 

For sufficiently large values of n and hence 
■y/— 6^ — 26 log A, we the second term on the right 
hand side can be ignored to obtain 



A(6,n)2 



e e 



b^-(logn)V46 



(29) 



and 



(35) 



where we have kept the root that is positive for all values 
of c. From Eqns. lT4 l l33l and 135 ] we can see that 



where 



(36) 



(37) 



P,{b,x) = e-'' e-('°s«)V4b^„_ (30) 

From the properties of Gaussian integrals it follows that 
if large, the above integral gives 



26e ^x ^-(iogx)V4b^ 
logx -26 



(31) 



Finally, from Eqns. lT2 l l33l and 135 ] we obtain the finite- 
entanglement scaling of the entropy, which is 



5* 



(v/f 



logx- 



(38) 



